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We recently proposed a simple dilaton-derived quintessence model in which the scalar field 
was non-minimally coupled to cold dark matter, but not to 'visible' matter. Such couplings can 
be attributed to the dilaton in the low energy limit of string theory, beyond tree level. In this 
paper we discuss the implications of such a model on structure formation, looking at its impact 
on matter perturbations and CMB anisotropics. We find that the model only deviates from 
ACDM and minimally coupled theories at late times, and is well fitted to current observational 
data. The signature left by the coupling, when it breaks degeneracy at late times, presents a 
valuable opportunity to constrain non-minimal couplings given the wealth of new observational 
data promised in the near future. 
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I. INTRODUCTION 

There is recent evidence ( [0- ||) that the Uni- 
verse's expansion is accelerating. If this is so, it would 
have fundamental cosmological implications, for pro- 
gressing the dark matter problem and reconciling a 
high Hubble Constant, h ~ 0.65, with an old Uni- 
verse to > HGyr. To explain such an acceleration, 
the Universe would have to have a matter component, 
additional to ordinary matter and radiation, since the 
latter two have equations of state that are unable to 
generate the required kinematics. In line with cur- 
rent observational constraints, the additional matter 
would have to have an equation of state p — wp with 

we (-i,-0.4)( [@- |). 

A pure cosmological constant cannot explain the ob- 
served acceleration without running into fine tuning 
problems; one would need A ~ 10~ 122 c 3 /(hG), sev- 
eral hundreds of orders of magnitude lower than one 
would expect from a vacuum energy originating at the 
Planck time j7|. This has lead to a wealth of propos- 
als using a scalar "quintessence" field, minimally cou- 
pled to matter through gravity, which can be cajoled 
into acting as an effective cosmological constant in the 
presence of a suitable potential. Models of particular 
interest use "tracker" potentials (e.g. [||- |l3|]) which 
allow the scalar field to produce the required dynam- 
ics without dependence on initial conditions, but these 
still require small-scale parameters. More recently, a 
model was proposed ( |Q, |l5|) with a potential who's 
parameters were, a more physically agreeable, Planck 
scale. Explaining why the acceleration has only arisen 
recently, however, still requires some degree of fine- 
tuning in the model parameters, if not in the initial 
conditions, in order to confine acceleration to the cur- 
rent epoch iflql . A more practical explanation for the 
coincidental acceleration nowadays is that we are in 
close proximity to the cosmological transition from 
radiation to dust domination. Armendariz-Picon , 
Mukhanov and Steinhardt fllTf utilised this proxim- 



ity to drive the dynamics of their K-essence model 
although the Lagrangian used is somewhat complex, 
consisting of a series of non-linear kinetic terms. 

In a recent paper | fl8[ , we proposed a simpler model 
which harnesses the dynamical shift in the radiation- 
dust transition using a non-minimally coupled scalar 
field. We showed that a coupling of this form can 
use the transition to dust domination to push a 
quintessence field off scaling behaviour, and produce 
acceleration in the background nowadays. 

In this paper we consider the impact of such a non- 
minimal coupling on the evolution of perturbations to 
the background and the subsequent implications for 
both structure formation and the Cosmic Microwave 
Background (CMB). 

We start by giving an overview of the coupled 
quintessence model, and then go on to discuss the im- 
plications of coupling for perturbation evolution and 
structure formation. 



II. COUPLED QUINTESSENCE MODEL 

Non-minimal theories are commonly expressed in 
one of two frames. In one, the problem is posed in the 
Jordan frame and the scalar field is directly coupled to 
curvature, in the form f((j>)lZ, and produces a depar- 
ture from Einstein's gravity, as is seen in Brans-Dicke 
theories . This effect was used by [^0| to force the 
quintessence field out of scaling behaviour, necessary 
to give accelerated dynamics, however this "R-boost" 
occurs early in the radiation epoch and cannot ex- 
plain acceleration today. In the second, the Einstein 
frame is used and the scalar field instead couples to 
terms in the matter Lagrangian resulting in dynami- 
cal, field-dependent, masses and polarisations. These 
two groups are interrelated through conformal trans- 
formation of the metric; any theory in one frame can 
be rephrased in the other. However, usually a sim- 
ple function in one frame is mapped into a compli- 
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cated function in the other. Such couplings are heav- 
ily constrained when applied to the visible matter in 
the Universe, whether to photons ]2~l] ], or to what is 
usually called baryons (2^]. However, it could be that 
the dilaton coupled differently to visible matter and 
to the dark matter of the Universe. This hypothesis 
was suggested in pj|, and allows for large couplings 
to be consistent with observations (see also |24| , |2j| ) . 
We consider a scenario in which such a case exists. 

We choose g^ v to have convention (+---) in a 
flat FRW background. All quantities are expressed in 
units with M P = (8ttG n )~ 1/2 = 1 where M P is the 
Planck mass and Gjv is the Newtonian Gravitational 
constant. We consider a Lagrangian of the form: 



£ 



R 
~2 



V(4>)+c v + f (<!>)£, 



(i) 



in which £y is the Lagrangian of "visible matter" 
(baryons, photons, and also baryonic and neutrino 
dark matter), and £ c the Lagrangian of a domi- 
nant non-baryonic form of cold dark matter. We take 
V(cf>) — Voe~ x ^ the standard quintessence potential, 
which drives scaling behaviour when the coupling is 
minimal ( Q, Q). 

The coupling investigated is of the form f((f>) = 
1 + a(<f) — </>o)' 3 - Couplings of this form could arise 
as generalisations of an effective action for massless 
modes of a dilaton |2^] after performing a conformal 
transformation from the string frame into the Einstein 
frame, a and (3 are parameters reflecting the shape of 
the minimum being approached by the coupling [ErJ . 



A. Background Evolution 

Here we discuss the background equations in the 
conformal FRW metric, which are derived in appendix 
VI A . The field equations are obtained by varying the 
action with respect to the metric and the scalar field: 



Gr = Tjp+T$ + f{4>)T$ 



(2) 
(3) 



where G Mt/ is the Einstein's tensor and the various 
T M „ are stress-energy tensors. Heuristically, we may 



interpret the new term driving <p 
an effective potential V e ff = V 
identity (V^G^ = 0) leads to : 

VjjT ^(v) = Q 



as a contribution to 
- f(<p)£ c . Bianchi's 



(4) 
(5) 



These are to be contrasted with Amendola's coupled 
quintessence p7j ] (for which the interaction term is 
proportional to T). 
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FIG. 1. The evolution of S}^ 
with A = 8, /3 = 8,a = 50 
Qb = 0.053, h = 0.65). An early period of scaling is broken 
near the transition from radiation to matter, first with a 
period of kination, then inflation. At late times the uni- 
verse returns to a matter dominated scaling solution. 



Evaluating the components of the field equations, 
with scale factor a, we find Friedmann equations: 
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Pb + Pi + /W)Pc + ™ +V(tf>) (6) 



( Pc + £ C ) = Q (7) 

(8) 
(9) 

2 = -}'{<f)Pca 2 (10) 



where dots represent derivatives with respect to con- 
formal time, and the prime (') indicates differentiation 
with respect to <f>. 

One notices in equation (^) that the evolution of 
the background coupled dark matter is unaffected by 
the coupling. This simply arises because we are cou- 
pling to pressureless matter for which £ c = — p c ; if we 
had instead coupled to radiation we would find £=0 
and the coupling would have altered the background 
evolution (as discussed in appendix VIA ). However 
as will be discussed later, observations measure the 
coupled energy density f{4>)p c not simply p c so that 
the magnitude of the observed matter is affected by 
the coupling through (110). 

Fig.l. shows the evolution of and overall equa- 
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FIG. 2. The amount of accelerated expansion produced 
with the model for various values of /3, measured in the 
number of e-foldings N e = af/ai. ai and a / are the 
expansion scales when inflation begins and ends, respec- 
tively, (ai < ao < a/, where ao = 1 is the expansion scale 
nowadays). We have given a = 25 as an example with 
4>o = 260/A. 
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FIG. 3. The amount of accelerated expansion produced 
with the model as a is varied, measured again in the num- 
ber of e-foldings iV e .We have given (3 — 6 as an example 
with 0o = 260/A. 



our discussion below. 



tion of state w to t = Ptot/ptot for one model scenario. 
One can see that deep in the radiation epoch the cou- 
pling has a negligible effect on the overall dynamics 
and the scalar field's energy density scales with that 
of the dominant radiation, as in the minimally-coupled 
case. As the transition from radiation to matter domi- 
nation is approached the coupling becomes important 
and the dynamics are driven away from scaling be- 
haviour. The driving term on the right hand side of 
equation (|l0|) first, transiently, drives the field to ki- 
nate, suppressing the evolution of the scalar field and 
Vttj, ~ , then it re-emerges into inflationary behaviour 
to provide the accelerated expansion we observe today. 
The model requires that (3 be even and that the value 
of (f>o is of the order of magnitude of the scalar field 
today. However given these constraints, the model 
provides acceleration for a wide range of parameters 
as shown in the parameter space plots for the non- 
minimally coupled model in Figs. 2 and 3. 

In minimally coupled models with exponential po- 
tentials, the value of the parameter A is limited by 
BBN constraints Q to be A > 8 however the NMC 
model avoids this constraint through the suppression 
of Sl^ at nucleosynthesis, irrespective of A's value. Pa- 
rameter constraints for the non-minimal case can only 
therefore come from CMB and matter power spectrum 
predictions discussed below. In order to compare the 
non-minimal models with analogous minimally cou- 
pled ones however, we consider cases with A = 8 in 



B. Observational implications 

It is pertinent to consider whether the effect of the 
non-minimal coupling on the background at late times 
could be seen in current observations, i.e. when look- 
ing at the predicted apparent magnitude versus red- 
shift relation at z < 2. 

The apparent bolumetric magnitude is given by 

m(z) = M + 5logd L (z) + 25 (11) 

where M is the absolute bolumetric magnitude, and 
dh is the luminosity distance in Mpc 

"' = (1 + =>l3&) (12 » 

In Fig. 4 we plot the effective bolumetric magnitude 
from the B-band filter, m e J^ (z), from the Calan Tololo 
pj| and SCP Q surveys and predicted m(z) curves 
for the non- minimally coupled model in Fig.l and a 
comparative A CDM model. The effective magnitude 
is obtained from the apparent magnitude after tak- 
ing into account the lightcurve width-luminosity cor- 
rection, galactic extinction and the K-correction from 
the differences in the observed R-band and restframe 
B-band filters |(| . Within current observational error 
constraints, the non-minimally coupled model cannot 
be distinguished from the ACDM model. Recently 
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FIG. 4. Plot of the effective bolumetric magnitude 
for the Calan Tololo (open diamonds) and SCP data 
points (solid circles) against redshift. The curves corre- 
spond to two models considered in this paper A CDM 
model (solid line) with Q c = 0.347, Q b = 0.053 and 
SI a = 0.6 and a non-minimally coupled model with 
/(</>)n c = 0.347, fi 6 = 0.053 and = 0.6 with model 
parameters specified in Fig.l. 



proposed observational projects (see for example po| ) 
may offer future hope to discriminate between the ef- 
fect of quintessence models on the background evolu- 
tion. 

In the remainder of the paper, we consider 
an alternative approach to distinguishing between 
quintessence models, through their effect not on the 
background but on the perturbations about it. 



III. IMPLICATIONS FOR STRUCTURE 
FORMATION 

The addition of the scalar field has implications 
for structure formation both due to the addition to 
the homogeneous background energy density, and sec- 
ondly by the generation and evolution of scalar field 
perturbations. The additional background energy 
density shifts the equality redshift and alters the angu- 
lar distance to the last scattering surface. The scalar 
field also introduces extra terms in the perturbed Ein- 
stein equations and opens up the possibility of isocur- 
vature perturbations evolving. 

We study the impact of these effects by calculating 
the linear perturbation equations and specifying the 
initial conditions. These are then evolved from early 
on in the radiation epoch when the coupling is unim- 
portant through to nowadays. The matter and CMB 



power spectra are then calculated and compared with 
those obtained with minimally coupled models and 
observations. 



A. Linear perturbation evolution 

We follow the approach and notation of Ma and 
Bertschinger |3lJ extended by Ferreira and Joyce 
for minimally coupled scalar fields. A simplified model 
containing no baryons is used for the discussion, al- 
though a full theory containing baryons and relativis- 
tic neutrinos is used to obtain the CMB and matter 
power spectrum predictions presented. The essential 
results are presented here, while a full derivation of 



the equations can be found in the appendix VI B. 

Consider perturbations to a flat FRW metric in the 
synchronous gauge, with line element 

ds 2 = a(r) 2 {-dr 2 + (<% + h , , ul.r' d.r' } (13) 

We will only be concerned with the scalar modes of 
the perturbation, for which we can parameterise the 
metric perturbation as 



d i ke lk - x [k i k j h{k,T) + (hkj 



(14) 



where h is the trace of the metric perturbation. To 
obtain the linear perturbation evolution equations we 
consider the perturbed Einstein equations 

k 2 r) - ^Hh = 4vrGa 2 (5T (15) 

k 2 r] = 4:7rGa 2 ikiST° (16) 
h + 2Hh - 2k 2 r] = -8TrGa 2 6Ti (17) 

h + 6f? + 2H{h + Qri) - 2k 2 i] = 24vrGa 2 (fc l L - -£«)£!• (18) 

3 J 

where is the traceless shear. Writing the pertur- 
bations to energy densities, p, pressures, p, and the 
scalar field , in terms of a homogeneous background 
plus a perturbation, we have 



p(x,r) = p(t)(1 + S(x,t)) 
p(x, t) = p(r) + 5p(x, t) 
$(a;, r) = 0(r) + ip{x, t) 



(19) 
(20) 
(21) 



The only perturbation in to be affected by the 
coupling is 5T® , the other perturbations are the same 
as for a minimally coupled model, 

<5T ° = -p 7 <5 7 - (iff + fS c ) Pc - + <pV) (22) 



ihST° = \p 1 Q 1 + \^ 2 V 



ST! 



(23) 
(24) 
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where 9 is the velocity divergence. The evolution 
equations of the density perturbations for radiation 
and the dark matter component arc th e final require- 
ment. One finds, as is shown in VI B, that the cou- 



pling does not effect the first order equation for the 
matter perturbation so that, 



4 2 • 

3 7 3 
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(25) 
(26) 



The spare degree of freedom in the synchronous gauge 
allows us to choose the background, synchronous coor- 
dinates. As is conventional, we do this by constrain- 
ing the dark matter field such that 9 C = 0, which 
fixes S c = —\h. We are now able to write down the 
perturbation equations for the non-minimally coupled 
system. 

S c + HS C H — —S c 



the non-minimally coupled case early on in the ra- 
diation epoch, where the couplings effect is unimpor- 
tant. When the field is driven off tracking, close to 
the transition from matter to radiation, we cannot as- 
sume this, however. During the period when tracking 
is broken, the scalar field is suppressed and f2^ ~ 
(see fig.l). In general, the non-adiabatic pressure per- 
turbation 8p non - a d is given by 



8p 



non—ad 
P+P 



0(n*)(<y 7 + ^) 



(31) 



Therefore, since the quintessence contribution to the 
total energy density is highly suppressed, the isocur- 
vature contributions will continue to be small away 
from tracking behaviour, around the transition time. 
It is only at very late times, after last scattering, when 
is no longer small, that the isocurvature perturba- 
tions may start to grow. For the following discussion, 
therefore, we only consider adiabatic perturbations. 
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= -3irQ 7 <5 7 - 2ip(t) + {afV - 
$ + 2H£ + [k 2 + a 2 {V" + f'p c )]<p 
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(28) 
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The non-minimal coupling introduces extra terms 
into the equations for matter and scalar field pertur- 
bations, altering the mass terms and source terms , 
the latter shown on the right hand side of the equality 
for clarity. The coupling will only affect the radiation 
perturbations indirectly through the background bulk 
(via H) and through 6 C . 

Deep in the radiation epoch, the coupling to dark 
matter is unimportant. The adiabatic perturbation 
evolution closely follows the power-law solutions for 
the minimally coupled model with an exponential po- 
tential as discussed by Ferreira and Joyce |T^| . The 
growing modes of <5 7 , S c and (p evolve oc r 2 



<5 7 = 



<P 
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lo = n 

* 5A 



(30) 



where C is an arbitrary normalisation constant. 

It's only at very late times, z <^ 2, that the cou- 
pled matter establishes itself as the dominant effect on 
growth. This is when we would expect the coupling's 
signature to start to be seen. 

So far only pure curvature (adiabatic) perturbations 
have been considered, however isocurvature pertur- 
bations might also exist in quintessence models fj^] . 
For this non-minimal model we believe that their im- 
pact is negligibly small. Isocurvature perturbations 
are known to be negligible in minimally coupled track- 
ing quintessence models. This will also be so for 



B. Implications for matter perturbations 

An important consequence of non-minimal coupling 
is that, when considering the coupled matter, it is the 
coupled energy density, fp c , that should be interpreted 
as the matter density measured in observations, not 
p c ; an analogous case is non-minimally coupled grav- 
ity, f(4>)lZ, in which we consider the varying gravi- 
tational field strength as the observable and not con- 
stant Newtonian gravity, Gn- So we are interested in 
the effective dark matter density S c 



S(fPc) 
fPc 



Sc 



L 
f 



(32) 



An insightful way to look at the coupling's effect on 
perturbation growth is by looking at its effect on the 
dimensionless growth rate 



l eff 



Sc 
T — 

Sc 



(33) 



In Fig. 5 the growth rate for one scale, k = 
O.lMpc -1 , is shown for various models, in each case 
h=0.65, Sl fc =0.053. A non-minimally coupled (NMC) 
model with = 0.6 and A = 8(/3 = 8, cj) = 32.5), is 
compared with a ACDM model, n A = 0.6, a sCDM 
model O c = 0.947, and an analogous minimally cou- 
pled (MC) quintessence model using the potential de- 
veloped by Albrecht and Skordis |j V = V e- Xl >'{A+ 
(</>-<M S ) with A = 8 (A=0.0I, B=2, </> = 32.5). For 
z > 2 the growth rates for the scalar field models do 
not differ greatly from that in the ACDM model. 

The addition of a scalar field or cosmological con- 
stant, with Qq — 1 fixed, will act to reduce Q c and 
therefore the size of the mass term in equation (p7j). 
This is the main factor responsible for the suppres- 
sion of growth at later times, rather than the non- 
clumping behaviour of the scalar field commonly cited 
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FIG. 5. Time evolution of the effective growth rate for 
4 scenarios (all h=0.65, fif, = 0.053), 3 of which produce 
acceleration today: 1)ACDM ft A = 0.6 (full line), 2) NMC 
model = 0.6, A = 8, /3 = 8, (po = 32.5 (long dash) , 3) 
MC model = 0.6, A = 8, A = 0.01, B = 2, O = 32.5 
(short dash), and one which doesn't: 4)sCDM model 
f2 c =0.947 (dot-dash) 



as the cause. Sub-horizon scalar field perturbations 
have oscillatory time evolution with decaying ampli- 
tudes, their contribution to the evolution of matter 
perturbations therefore is small for the observation- 
ally interesting scales. For NMC models, the coupling 
suppresses Sl^ around z eq , making the scalar field con- 
tribution to 5 C growth negligible. Subsequently, the 
growth rate for NMC models is closer to that created 
by a cosmological constant than for the MC models. 

At late times however, for z < 2, the coupling and 
scalar field become important, and act to suppress the 
growth in S c to a far greater extent than A and MC 
models, offering a potential way to distinguish non- 
minimal from minimal the NMC model. 

The dampening effect can be also seen in the matter 
power spectrum P(k), 



P(jfc) = (\SJk)n = (100C)87r 6 h 6 k — 

k 



(34) 



where C is the normalisation factor from CMBFAST 
|3^| arising form the Bunn & White normalisation [Q 
at 1=10 multipole, k is in units of h/Mpc and ko = 
0.05Mpc _1 and n is the tilt, chosen here n=l for a 
scale invariant spectrum. 

In fig 6. the matter power spectra for the both the 
NMC and MC models mimic a ACDM model for scales 
k <~ 0.1 i.e. those modes having entered the horizon 
before and around equality. There is a slight suppres- 



k (h/Mpc) 

FIG. 6. Matter power spectrum for the 3 scenarios in 
Fig. 5 which produce acceleration today: ACDM (full), 
NMC (long dash) and MC (short dash), together with the 
de-correlated data points of Hamilton et al. Parameters 
in the 3 models are the same as in Fig. 5. 



sion but a bias factor could in theory resolve the dis- 
crepancy. Certainly all three models give reasonable 
predictions for matter fluctuations over a sphere of 
size 8/i- x Mpc with <r$ =0.89,0.91,1.13 for NMC, MC 
and ACDM models respectively, in comparison to the 
observed value <r s = 0.56ft" - 47 ~ 0.9 g§. 

For larger scales, however, the coupling does make a 
difference. In scales that have only entered the hori- 
zon in recent times, whilst the coupling is important, 
we see a distinctive reduction of power in comparison 
to the MC and ACDM models which tend to similar 
behaviour. Although the suppression clearly distin- 
guishes the coupled model, its profile is still consistent 
with current observational results ||(| . There may be 
an opportunity with the future SLOAN galaxy sur- 
vey results to constrain the power spectrum at these 
larger scales (to k ~ 0.01). 

Another potential impact of the late time impor- 
tance of the coupling is that it will affect small scale 
features at z ~ 2, observable potentially through 
future weak lensing (see e.g. |37j] and references 
therein) and damped Lyman a cloud measurements 
(see e.g j§). 



C. Impact on CMB anisotropies 

Introducing a scalar field can potentially have sev- 
eral effects on the CMB power spect rum. Firstly, as 
we have already mentioned in section HI A , the scalar 
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FIG. 7. CMB power spectra showing low 1 (plateau) be- 
haviour for the 3 scenarios in Fig. 5: ACDM (full), NMC 
(long dash) and MC (short dash) with COBE datapoints. 
Model parameters are the same as in Fig. 5. The 3 models 
evolve differently at late times producing slightly different 
ISW anisotropics shown in the plateau at low 1. However 
observations at this scale are dominated by cosmic vari- 
ance, so that the differences would not be observable. 



field gives rise to extra mass and source terms in the 
linear evolution equations for ip and 5 C . These then 
indirectly affect the radiation perturbations, altering 
the acoustic peak positions and heights at the time of 
last scattering (r; ss ). However, the scalar perturba- 
tions are effectively negligible around z eq , especially 
in the NMC scenario so this effect will be minimal. 

Secondly, the time varying Newtonian potential af- 
ter decoupling will be affected by the coupling, alter- 
ing the anisotropies produced at large angular scales 
(the Integrated Sachs- Wolfe effect). This can be seen 
in fig. 7 where the NMC model has a different profile 
at small 1 from the MC and ACDM models. However 
the effect is not large enough to be disentangled from 
the effect of cosmic variance. 

Thirdly, the inclusion of the scalar field alters the 
composition of the energy density, altering the angu- 
lar diameter distance of the acoustic horizon size at re- 
combination. This can be parameterised by the value 
of B 

B = nl /2 h dzlx^z^+^H (35) 

-'Zrec 

Altering the value of B shifts the positions of the 
peaks. The critical problem one confronts when try- 
ing to use CMB spectra to differentiate between mod- 
els is the degeneracy that exists between models with 
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FIG. 8. CMB power spectra showing acoustic peaks for 
the 3 scenarios in fig.5:ACDM (full), NMC (dash) and 
MC (short dash), together with the data from Boomerang 
(solid circles), Maxima (crosses) and DASI(open dia- 
monds). Model parameters are the same as in Fig. 5. 



identical fl c ,^b and B pS| ]. It has been shown that 
this degeneracy can be broken for scalar field mod- 
els in which a large fraction of the energy density 
at Tiss is from the scalar |^8[ ; the scalar field act- 
ing as an effective increase in the number of relativis- 
tic degrees of freedom. However for models in which 
A > 8 the degeneracy still exists in minimally coupled 
models. In Figs. 7 and 8, CMB spectra are plot- 
ted for the scenarios discussed in the previous section 
against COBE @, MAXIMA Q, Boomerang @ 
and DASI jl3| data. All the models discussed have 
23=1.77 and yet one can see that the degeneracy of 
the first peak is slightly broken, with the NMC model 
having Z pea fc=215 in comparison to 224 for both the 
MC and ACDM models. It is also interesting to note 
that the CMB spectra for coupled models with differ- 
ent A values are effectively degenerate in themselves, 
as shown in the figure. This implies that, although 
coupling itself may be distinctive, CMB spectra will 
not be able to isolate the parameter in the potential. 

The fourth possible effect is on the separation of the 
peaks. This has been proposed as a possible mech- 
anism with which to distinguish minimally coupled 
models [Q . They are not distinguishable from ACDM 
if fl^^Tiss) is small however, as mentioned above. But 
in the case of non-minimally coupled models the de- 
generacy in the second and third peaks is broken be- 
cause of the effect the coupling has on tq the con- 
formal time nowadays. This is of particular interest 
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given the expected improvements in peak definition 
(e.g H|, Hi), g|, H) . The separation of the peaks 
dl is given by 



51= IT 



TQ — Tls 



(36) 



where r s is the sound horizon and c s , the baryon speed 
of sound, both of which can be assumed effectively 
constant across the models. The NMC model has 
to=12,530 in comparison to r =13,077 for the ACDM 
model. This reduces the separation slightly breaking 
the degeneracy, as shown; the separation of the first 
and second peaks in the NMC model is dl=309 in 
comparison to 327 ACDM scenario. 

Although distinguishable from the cosmological 
constant spectrum, the difference is still too small to 
be resolved with current observational data, including 
the most recent Boomerang [Q and DASI Q] data, 
showing highly improved definition in the second and 
third peaks. However with a number of observational 
projects continuing to focus attention on resolving the 
higher peaks, the breaking of degeneracy may offer a 
way to constrain non-minimally coupled models. In 
Fig. 9 we plot the residual differences between the A 
CDM and NMC Ci spectra in Fig. 8 when compared 
with estimated MAP errors. The parameters used to 



estimate the MAP errors are shown in appendix VI C 



The estimated errors are considerably smaller than 
these residual levels for I < 900 implying that we may 
be able to distinguish between these various models 
within the near future. 



IV. CONCLUSION 

We have examined the impact a scalar field, non- 
minimally coupled to cold dark matter will have on 
the evolution of matter and radiation perturbations. 
We considered firstly its impact on the linear evo- 
lution equations and found that even though it did 
introduce new terms these were effectively negligible 
for all but very late times. The impact of this late 
time behaviour was then considered for matter per- 
turbations where it was seen to create a suppression of 
growth at large scales. The coupling was also found to 
break the degeneracy usually seen in the CMB spec- 
trum, slightly shifting the position of the first peak 
and reducing the separation between adjacent peaks. 
These two distinctive 'signatures' of the coupled dark 
energy model are not resolvable with current observa- 
tions. However projects currently underway look to 
mapping both the matter power spectrum and CMB 
peaks with much improved accuracy. These may of- 
fer an opportunity to eventually distinguish between 
ACDM, minimally coupled and non-minimal coupled 
quintessence models in the near future. 
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FIG. 9. Residual deviation of the the NMC model's 
temperature fluctuations from those of the ACDM model, 
both shown in Fig. 8., with estimated MAP error bars. 



VI. APPENDIX 

A. Bianchi's Identity 

In this section we derive the equations of motion 
for the coupled dark matter and scalar fields explicitly 
from the action. 

S = [ ^d 4 x(-^ + l< 



^d i x{-- + -4>^ -v{4>) + c v + f{4>)Cc 

(37) 

Bianchi's identity reflects the symmetry of the Rie- 
mann tensor, the Einstein tensor being convariantly 
conserved, 

_ -g^TL).^ = (T^ w + T^ v{v) + f{(j))T^ {c) ), u = 

(38) 

Since visible matter is minimally coupled in the model 
we can immediately separate it out, 



T! *v{V) = o 



(39) 



Using the explicit definition of the energy momentum 
tensor in terms of their Lagrangian, 
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rj-ijXV 



-9 



(40) 



so the scalar field component of the Bianchi identity 
in terms of 6 and its derivatives is 



i™ = <K {K + KA a + V {<!>)) 



(41) 



The Euler-Langrange equation, which is just the Klein 
Gordon equation, allows us to simplify this further 



Combining these two expressions we obtain 
For the coupled matter then, 

Combining the results in equations (|39| ) and 
Bianchi's identity in (Bq) is given by 



^j4>;v {f-c9 



(42) 
(43) 

(44) 
©, 

(45) 



We can obtain an expression for the Lagrangian for 
perfect fluid by considering it to be a gas of particles 
with masses m a and paths x l a f4(| 

C a (x) = -m a S(x - x a (t)) {~g^x v y (46) 
By noting that the length of the 4-velocity 

(-g^iirijj:) 2 = ds/dX equals 1 for dust, this ex- 
pression simplifies greatly. Averaging over particles 
in the gas rest frame we find C c = —n(m/u°) where n 
is the particle number density and m m = dx^/dX is the 
4-velocity. For pressureless particles — {1,0,0,0} 
therefore L c = —p c . 

We can also obtain an expression for the stress- 
energy tensor from the Lagrangian in an analogous 
way . Using the relationship between T^ v and L in 
(40) the stress-energy for a particle is given by 



rr=m/ (a ^ (t) W u ° 



(47) 



On averaging over the particles the energy density 
is p = (T 00 ) = n(mu°) and the pressure in the x- 
direction is given by p = (T 11 ) = n(mu° (v 1 ) 2 ) where 
v i = dx l jdt = u l /u°, for dust therefore p c = 0. In 
addition, in the rest frame there is zero streaming ve- 
locity so that (T 0i ) = n(muV) = 0. 

In this paper we assume that the coupled matter is 
comprised of cold pressureless dust particles. Putting 
the expression for C c into (|45|) , the background evo- 
lution equation for non-minimally coupled matter is 
identical to that in the minimally coupled case. 



B. Linear perturbation equations 

We here derive the linear perturbation equations for 
the coupled cdm and scalar fields in detail. We assume 
the notation of Ma and Bertschinger |3l| and results 
of Ferreira and Joyce E3] . Consider perturbations to 
a flat FRW metric in the synchronous gauge, with line 
element 

ds 2 = a(r) 2 {-dr 2 + (% + //,, .d.r'il.r' } 

By considering the perturbed Euler-Lagrange equa- 
tion for the scalar field we can obtain the linear evo- 
lution equation for (p. 

< + + r%v a + v"y = f"v£ c + f'sc c (49) 

Perturbing the particle Lagrangian in ( fl6"l ) 

8C a = ~Sm a 8(x - x a (t)) {-g^x v y (50) 

and averaging over all dust particles, we have SC = 
—n(6m/u°) = —Sp c . The non-minimally coupled 
scalar perturbation equation becomes 

^ + + 2% + V"<p - = -p c (/'V + AO (51) 

The perturbed Einstein equations can be used to ob- 
tain an expression for Sp 

(p + P)9 + 3H(5p + 5P) + ^{p + P)h + 6p = (52) 

We are interested in the interacting dark matter and 
scalar field, which are not separable in (|5^ ) 

fp c 9 c + 3HS(fp c ) + \{fp c )h + -^(S(fPc)) 



a 2 ~^ a 2 



V'ip + V'\ 



(53) 



Using the equations of motion for the background and 
perturbed scalar field we find that (^3|) simplifies sub- 
stantially. Interestingly, we find that the coupling 
does not affect the first order dark matter perturba- 
tion equation 



5c = -0c - \h 



(54) 



(C) = Pc + 3H Pc = 



(48) 



With the residual degree of freedom in the syn- 
chronous gauge we are free to fix one additional pa- 
rameter, by convention we set 8 C = so that 5 C = ^h. 
Ignoring baryons, for simplicity, the second order per- 
turbation equation becomes, 

S c + HS e + ^6 C = 

-3H 2 n 7 S 7 - 2# + (oV - 3H ^' )tp (55) 
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C. MAP error bar estimation 

The standard error on the estimate of Ci, ACi for 
an experiment with N frequency channels (denoted 
subscript n) , each respectively with angular resolution 
0n,fwhm(o,rcmin) = n (rad)/6O x 7r/180 and sensitiv- 
ity a n per resolution element, scanning a fraction f s k y 
of the sky in bins of I size Al is given by |l7| 



AC; 



0.5 



[Ci+u-iBf*] 



{2l + l)f sky Al 



U).„ 



(cr n O n )~ 



B 



nl 



-l{l+l)/l 2 B 



(56) 



where we have assumed that the experimental beam is 
approximately Gaussian filtering scale l s = y/8ln29~ 1 . 
We have assumed a useful sky fraction f s k y — 0.65 
and Al = 50, Table 1 shows the remaining parameters 
used, taken from pfifl . 

TABLE I. MAP CMB Experimental Specifications 



v (GhZ) 


@n,fwhm 


10 6 cr 


N ch 


40 


28' 


8.2 


4 


60 


21' 


11.0 


4 


90 


13' 


18.3 
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The combined errors in Fig. 9 are 
therefore A[Q{ACDM) ~ Ci{NMC)]=ACi{KCDM) 
+ACi{NMC). 
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